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ABSTRACT 

A stabilization problem for Burgers’ equation is considered. Using linearization, vari- 
ous controllers are constructed which minimize certain weighted energy functionals. These 
controllers produce the desired degree of stability for the closed-loop nonlinear system. A 
numerical scheme for computing the feedback gain functional is developed and several nu- 
merical experiments are performed to show the theoretical results. 
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1. Introduction 


Burgers’ equation is a one-dimensional simple model for convection-diffusion phenomena 
such as shock waves, supersonic flow about airfoils, traffic flows, acoustic transmission, etc. 
For high Reynolds number, the solution of Burgers’ equation produces steep gradients due 
to the nonlinear nature of the convection. 

In this paper, we consider a stabilization problem for Burgers’ equation. We employ 
“a-shifted” linear feedback control laws to obtain a desired degree of stability, on a certain 
energy space, for the closed-loop nonlinear system. 

Burgers’ equation 

^ ^ q2 

(1.1) + Z(t,x) — z{t,x) = 6 — 2(t,a0 

was introduced by Burgers [5,6,7] as a simple model for turbulence, where e > 0 is a viscosity 
coefficient. Since then, many researchers have considered the conservation law 

3 3 

(1.2) —z{t,x) + z(t,x)—z{t,x) = 0 
and the “viscosity solution” 

(1.3) z{t,x) = lim z c (t,x), 

v ' ej.0 


where z t (t,x) satisfies equation (1.1), see [9,16,18,22,24,25,26]. 

Oleinik [26] proved that for any T°°-initial data, there is a unique viscosity solution for 
equation (1.2) and the solution satisfies the “entropy condition” 


(1.4) 


z(t, x + a) — z(t, x) E 
a t 


for all t > 0, a > 0, -oo < x < oo and for some constant E > 0. A complete discussion of 
these results may be found in [32]. 


Almost no results exist for the control problem associated with Burgers’ equation. Chen, 
Wang and Weerakoon [8,35] considered an optimal control problem for equation (1.2) with 
-oo < x < oo. The problem was to select an initial function to minimize a specific cost 
functional J. They obtained sufficient conditions for the differentiability of J with respect to 
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the initial function and explicit expression of the entropy solution of (1.2) in terms of initial 
data. 

In this paper we consider a control problem for Burgers’ equation (1.1) defined on a 
finite interval. Specifically, we will find several feedback laws stabilizing the nonlinear system 
(1.1) with a certain exponential decay rate. The feedback laws will be obtained from the 
linearized equation. Curtain [11] has considered a stabilization problem for certain semilinear 
evolution equations. Using Kielhofer’s stability results for semilinear evolution equations [20], 
she showed that, under certain conditions, there exists a finite dimensional compensator 
which produces a stable closed-loop system. These finite dimensional compensators are also 
obtained from the linearized control system. Applying her results to Burgers’ equation (1.1) 
with, for example, Dirichlet boundary conditions, one can obtain the stabilizability results 
of the closed- loop system which are similar to ours. However, in [11], there is a restriction on 
the action of the output operators. The domain of the output operator was required to be 
a certain subspace of L 2 which contains the Sobolev space Hq. In this paper, we investigate 
optimal feedback laws, instead of finite dimensional compensators, in the sense that they 
minimize certain energy functionals. 

Well-posedness and stability results for the open-loop system are obtained in Section 2. 
In Section 3, a “shifted” linear control problem, (LQR) 0 , is introduced. Under appropriate 
selection of input and output operators, (LQR) * is stabilizable and detectable. The feedback 
control law obtained from (LQR)q produces the desired degree of stability for the closed- 
loop nonlinear system, (Theorem 3.10). In Section 4, a numerical scheme for computing the 
“feedback functional gains” is developed and several numerical experiments are performed. 

We shall use standard notation. If (X, || • || x ) and (Y, || • ||y) are normed linear spaces, 
then C(X, V ) will denote the space of all bounded linear operators from X to Y and for any 
A G C(X, Y), ||A|| or ||A||£(x,k) will denote the operator norm on the space C(X : Y). In 
the event that X = Y we denote C{X,Y) by C{X). From time to time we will use || • || 
without any subindex for vector or operator norm. In all such cases the appropriate index 
for || • || will be understood from the context. For a Hilbert space X, we denote the inner 
product on X x X by < •,• >*. Given a linear operator A from X into itself, we denote 
its domain, spectrum, resolvent and adjoint by V(A), cr(A), p(A) and A *, respectively. For 
real numbers a, b with a < b, L p (a,b‘,X), 1 < p < oo, will be the space of all Lebesgue 
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measurable functions / from (a,b) to X such that ||/||z,p(a,&) = (/„ \f( x )\ p dx)r < oo. The 
spaces H k (a,b ) and H^a^b) are the standard Sobolev spaces defined by H k (a,b ) = {/ £ 
L 2 (a,b)\f^ e L 2 (a,b),j = 0, 1, — ,k) and H£(a,b) = {/ £ H k (a,b)\f^(a) = / (j) (6) = 
0, j = 0, 1, • ■ • , k — 1}, respectively. The dual space H~ k (a,b) of H£(a, b) is the space of all 
continuous linear functionals on H^a^b) represented by the inner product < >L2( a ,6)- 

Finally, we present a physical example that may be found in most standard references to 
motivate the control problem for Burgers’ equation. Other examples involving traffic flows, 
supersonic ffow about airfoils, acoustic transmission and turbulence in hydrodynamic flows 
can be found in [12] and the references given there. The following example is taken almost 
directly from [12]. 

Example (Shock Waves) 

An impulsively-started piston moving at a constant velocity into a tube containing a 
compressible fluid initially at rest creates compression waves. The compression waves even- 
tually coalesce, due to the nonlinear nature of the convection, to form a single shock wave. 
The one- dimensional unsteady motion of the fluid is governed by the continuity equation 

d d d 

(1.5) dt P ^ t,xS> + + v { t i x )jfaP( t i x ) = 0 

and the x-momentum equation 

d d d ^2 

C 1 - 6 ) di v ( t ' x ) + x ) + X ))M*> x ) = x )’ 

where p is the density, v is the velocity, p is the pressure and 6 is the “diffusivity of sound”. 
It is convenient to replace the density by the sound speed, a = a(t,x ) by a(t,x)/a 0 = 
(p(t,x)/ po) 2 ^- , where 7 > 1 is the specific heats ratio and the subscript 0 refers to the 
undisturbed values [23]. Equations (1.5) and (1.6) become 

(1-7) jU(t,x) + v(t,x)-^-a(t,x) + ?-y^a(t,x)-^v(t,x) = 0 

and 

d d 2 c) ^2 

( L8 ) di V ^' X ^ + + ’^T\ a ( t ' x )fr c a ( t ' x ) = ^^ 2 U ^’ X )’ 
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where 8 is a function of the undisturbed (to the right of the shock) viscosity, density, specific 
heat and thermal conductivity of the medium. Equations (1.7) and (1.8) can be simplified 
by introducing the Riemann invariants, 


(1.9) 


to give 


( 1 . 10 ) 


and 

(i.ii) 


r(t,x) 


a(t,x) v(t,x) 

7^T + ~2~’ 


s(f,:r) 


a(t,x) v(t,x ) 

7 - 1 2 


d . .. 5 . . 8 d 2 

— r(t,x) + (a(«,z) + v(t,x))—r(t,x) = -^( r (t,z) - s(t,x)) 

$ d & ^ 

— s(t,z) - (a(t,x) - v(t,x))—s(t,x) = -^2(s(t,z) - r(<,*)). 


Consider the propagation of a disturbance into an initially undisturbed region, s = Sq 
where sq = Then the problem is governed by equation (1.10). But from equation (1.9), 


( 1 . 12 ) 


a{t,x) + v{t,x) = - ^ ^ r(/.x) + — 




thus equation (1.10) becomes 

(1.13) j t r(t,x) + (^(M) + = ~r((,x). 

As the final step we introduce the change of variables 

(1.14) z(t,x) = -^Y~(r(t,x)-r 0 ), f = x - a 0 t 
to give Burgers’ equation 

(!- 15 ) £*(*>0 + = 

From equations (1.12) and (1.14) we have z(t,£) — (a(t,£) + u(f,£)} — { v 0 -f a 0 }, where 
u 0 = 0, i.e., u(£,£) is the excess wavelet velocity (the difference between propagation speeds 
of disturbance in stagnation and nonstagnation conditions). The coordinate £ is measured 
relative to a frame of reference moving with the undisturbed speed of sound a 0 . 
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2. Well-posedness and stability of Burgers’ equation 

In this section, we consider well-posedness and stability properties of the solution for 
Burgers’ equation with Dirichlet boundary condition. These results will be needed in the 
analyses of our control problems in the next section. We first consider an abstract version 
of this problem and then specialize to Burgers’ equation. 

Consider an initial value problem 

(2.1) ^z(t) = Az(t) + f(t,z(t)), z(t 0 ) = z 0 , (t>t 0 ), 

on a Hilbert space X , where A is the infinitesimal generator of an analytic semigroup S(t ) 
satisfying ||5(<)||£(x) < Me“‘, t > t 0 , for some constants M = M(u) > 1 and to > 0. Since 
S(i) is analytic, the fractional powers of Ai = -A + al are well-defined for any a > to [27, 
Chapters 1,2]. Since 0 £ p{A\), the resolvent of A u is invertible for all 0 < \i < 1. 
Therefore, the graph norm 1 1 a: 1 1 4- ||^4 j 2|| on the domain V(Ai) of Af is equivalent to the 
norm ||z|| M = ||.45*z||. We denote the Hilbert space V(A^) with the norm ||z|| + \\A^z\\ or 
|Kz||by(X M ,||.||,). 

We shall make the following assumption [17,27]. 

Hypothesis (F): Let U be an open subset of [t 0 ,oo) x X The function / : U — > X 
satisfies the hypothesis (F) if for every ( t,z ) £ U there is a neighborhood V C U and 
constants L>0, O<0<1 such that 

(2.2) ||/(<i,^i)-/(< 2 ,^)IU < L(\t i-t 2 f + \\z x -ar 2 |y 

for all (<,-, Zi) € V, i — 1,2, i.e., / is locally Holder continuous in t , locally Lipschitzian in z , 
on U. 

Now we are ready to state the local existence theorem for the solution of equation (2.1). 
The following result appears as Theorem 3.3.3 in [17]. 

Theorem 2.1. Let A be as before and / satisfy hypothesis (F). Then for any (to,z 0 ) £ 
U C R + x Xp, there exists T - T(t 0 ,z 0 ) > 0 such that equation (2.1) has a unique (strong) 
solution z(t) on [t 0 ,to + T) with initial value z(t 0 ) = z 0 - 

Now, consider Burgers’ equation, with Dirichlet boundary condition, on a finite interval 
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[0, £} given by 


^ q2 ^ 

*^) ^(^7 ®) ^(^7 *^)? 0 < X < f , i >0 

z(<,0) = z(t,^) = 0 

z(0,x) = 2o(x), 

where e = ^ > 0 and Re is the Reynolds number. In order to place the system (2.3) into 
a semigroup framework let z(t)(-) = z(t,-), z 0 (-) = z(0 , -) and H = L 2 (0J). Define the 
operator A c by 

(2.4) A c <f> = e (/>" 

for all <f> G 'D(A C ) = H 2 (0,£) D Hq(0,£). The system (2.3) can now be written as the initial 
value problem 

(2.5) j t z{t) = A t z(t) + f(t, z(t )), z( 0) = z 0 , ( t > 0), 

on the space H, where f{t,z) = —zz‘ is defined on the space Hq(0,£). It is well-known 
[17,27,34] that A t generates an analytic semigroup S(t ) on H. 

We summarize the basic properties of the infinitesimal generator A ( and its semigroup 
S(t), t > 0, in the following remark. 

Remark 2.2. (i) The spectrum <r(A t ) of A e consists of all eigenvalues A„ = -en 2 7r 2 /^ 2 , 
n = 1,2, • • • , and for each eigenvalue A n the corresponding eigenfunction <f> n is given by 

(2.6) <f>n(x) = V2 sin ~j-x, 0 < x < £. 

(ii) The operator A t is self-adjoint, i.e., A c = /!*, and the semigroup S(t) can be repre- 
sented by the following formula 

OO 

(2-7) S(t)z = £ < z,<j> n > <j> n 

n= 1 

for all z £ H, where <f> n 's are defined by equation (2.6). Moreover, from equation (2.7), it is 
easy to see that S(t ) has the stability property 

(2-8) ||SW|| £( h> < e-l"’" 1 )', t > 0. 
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A simple application of Schwartz inequality gives the following first Poincare inequality 
[36, p. 116]. 

Lemma 2.3. For any z € Hl(0,C), 

(2.9) 11*11 H < *ll*'lltf» 

where H = L 2 { 0,i). 

Remark 2.4. (i) The above lemma gives an equivalent norm ||^|[//i = \\ z '\\l 2 on the space 

^(0,0- 

(ii) It is well-known that — Hq(0,£) [17, p. 29], [21, p. 326]. 

Lemma 2.5. For any 2 G Hj(0,£) = V((-A e )s), the following inequalities hold. 

(2.10) ||5(0*lla8 <^11*11^, (*>0) 

(2.11) ||S(< Wk <(^i + ^f5)e-’ t |W| i =(o,,). (*>"). 


where 7 = e7r 2 /^ 2 . 

Proof: For any z G Hq( 0,£), we know that (-A e )iS(t)z = S{t)(— A c )*z. Hence, by 
Remark 2.2, 

||S(f)z|| H i = l|S(t)*lk* + ||(-A e )i5'(<)*||L» < 11^(011 (11*11 + ||(-^e)MI) < e 7< ll*ll«J- 

The inequality (2.11) follows from Remark 2.2, Lemma 2.3 and the estimate 

OO 

l|s(i)z|| H . < (i +/) ||(-4)ls(«)z||w = (i + 0 ll(-A) 4 Y, eK ‘ < 2 ^» > *»H 

n=l 

co 

= (i + i )\\^2-f eXnt <z ^ n> ^ cos ~f x \\ 

n = 1 

< (1 +£) (sup{^e (An+7)t : n = 1, 2 , . . . })e" 7t ||z||// 


and 


su P{-^“ e 




n 



0<*< 5& 


t > 


e 

2 CK 2 ’ 


where A n = — en 2 7r 2 /^ 2 and <f> n ( x ) — \/2sin ^ x, n 1,2,***. 
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Remark 2.6. The inequality (2.11) holds for every z £ H = L 2 ( 0,£), since the semigroup 
S(t ) is analytic. 

Now we have the well-posedness and stability properties of Burgers’ equation (2.3) on 
the space H£(0,l). The following theorem follows from an application of Theorem 5 1 1 in 

[171- 

Theorem 2.7. For any given /? > 0, 0 < /? < 7 = en 2 /£ 2 , there is a p = p(£,e,0) > 0 
such that for any initial data z 0 6 0,1), with ||z 0 ||„. < §, there is a unique solution 

z{i) = z(t, 0; Zq) E H£(0,£) of equation (2.3). Moreover, the solution satisfies the inequality 

(2-12) lk(< ? 0; zo)\\ H \ < 2e — ^ll^ollz/i (t > 0) 

and p = p(t , e, 0) > 0 can be chosen to satisfy 

(2-13) 0 < p < — ZR 

V^TT (1 + £){l\J 7 — /? + y/27re ) 

Proof: Note that 0 is an equilibrium point for the system (2.5). Since £>((— yj e )|) = Hq 
if the nonlinear term f(z) = -**' satisfies the hypothesis (F) with index p = I, then, by 
Theorem 2.1, we have a unique local solution z(t,0;z o ) on the space H r 0 x . It is easy to see 

that | \f(zi) - f(z 2 )\\ L 7 < (lki||// o i + INi - for all z u z 2 E tfj, uniformly in 

i > 0. Hence, / satisfies the hypothesis (F). 

For the global existence and uniqueness of the solution z(t,0;z o ) E let z 0 be any 

initial data in H l 0 with ||z 0 ||//i < f, where p = p(£,e,/3) satisfies the condition (2.13). It 
follows that 


(2.14) 



V5 


e -bt-P)s ds + 


(! + ^)?r 

% - 0 ) } 



? 


s * nce Jo ds y^Zp r(i) sinf — w ^ ere r is the Gamma function and T(|) = 

Let 1 1 1 1 // * < 2 ‘ ^hen, by the local existence property, there is a unique solution 
z(t,0;zo) £ Hi satisfying the inequality \\z(t,0; z 0 )\\ H i < p on an interval [0,< x ) for some 
t\ > 0, where t x is chosen as large as possible. We will show that t x must be infinity. Suppose 
that t x is finite. Then we must have ||ar(i 1 )|| H , > p. Note that, on the interval [0,t x ), 

^ 2 ’ 15 ^ ll/(*(0)lli* = || - z{t)z'{t)\\ L 2 < (||a(0l|j/i) 2 < P 2 , 
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where ' = 4~. Lemma 2.5, Remark 2.6 together with inequalities (2.14), (2.15) yield 

ax 

rt\ 

ll^Olltfi = \\S{h)z 0 + / S{t 1 - s)f{z{s))ds ||„1 
0 Jo 

< ||5(<i>oIIh‘ + j 11^1 - *)/(*00)lk ds 

t , f r [ X - + (1+ -^ } ds 

Jo \/2e y/ti - 5 ^ 




i ±£ ri 
\/2e Jo V~s 


P 2 r l+< f°° 1 .-7* 

< - 4- p 2 r ' " 

”2 P 


e ds + 


(1 + ^)tt 

i 


roo 

/ e' 78 ds } 

Jo 




This is a contradiction. Therefore, the unique global solution z(t,0]z o ) exists. Moreover, 
from the above estimate, we know that if ||zo||tfi < f then ll 2 (*’ 0 o)||h> < P for a11 

t E [0, oo). 

Finally, we will derive the stability result (2.12). Let w(t) = sup{ ||z(s)||/ / ie /3 . 0 < 
s < t }. We then have 


Uz^lltfie* < + I \\S{t ~ s)f{z{s))\\ Hl ds) 

J 0 

< (e-|| 2 „||„, + l+i^L= + IW.)lli, * > 


i + i i , (i + 0* 


)e-^ t ^\\z(s)\\ Hl e^ds 


< IM„. + ^ ( jf * ) -(0 ^4 ^ - ) - w 

< IMk + H jf * + JfTfl } m(<) - IWI "" + 5“ (,) ■ 


Therefore, w(t) ^ 2||zo||//^ an< ^ ll z (0lltfd — || 2 o||//*- 

Remark 2.8. Rankin [29] considered well-posedness properties for a certain type of serru- 
linear evolution equations where the nonlinear terms are in divergence form. According to 
his results, we can see that equation (2.3) has a unique (strong) solution for initial data in 
L p (0,£), p > 4. To get this result, he used the analyticity of the semigroup S(t ) and the fact 
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that the differential operator £• on can be represented by 

(2.16) V~t±=(-A,)iB 

for some bounded operator B G C(H ), where H = L 2 ( 0,£). In general, it is not true that 
dx ~ ( ^t) 1 • This result could be used to analyze the stability property of the solution 
of Burgers’ equation with initial data in L p (0,(), p > 4. 

3. Linear control problem 

As we noted in Section 2, the open-loop (no control) solution of Burgers’ equation (2.3) 
decays exponentially in the topology of the energy space Hq( 0,^) (see Theorem 2.7). How- 
ever, the decay rate depends on the viscosity e > 0. We now explore the possibility of 
obtaining an exponential decay rate independent of viscosity by feedback laws. 

The basic model is governed by an abstract system of the form 

d 

( 3 J t z (t) = Mt) + Bu{t), z{ 0) = z 0 e H, 

y{t) = Cz{t), t > 0, 

where H, U and Y are Hilbert spaces, u(-) G T 2 (0, oo; {/), y (- ) G L 2 (0,oo;Y), and A is the 
infinitesimal generator of an analytic semigroup S(t ) on a Hilbert space H. Assume that 
BeC(U,H), C € C(H, Y) and A is self-adjoint with compact resolvent. Solutions of (3.1) 
are taken to be in mild form, i.e., 

(3-2) z{t) = S(t)zo + f S(t — s)Bu(s) ds , 

Jo 

(3-3) y(t) = CS(t)z 0 + C f S(t — s)Bu(s) ds } t > 0. 

Jo 

We first consider the performance index 

roo 

(3-4) J(u) = / { ||y(0lly + } dt, 

Jo 

where y(t) is given by equation (3.3), and R > 0. The linear quadratic regulator problem is 
(LQR) : Find u(-) G L 2 (0,oo;U ) minimizing the cost functional J given by equation 
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(3.4) subject to the system (3.2)-(3.3). 

For the existence of an admissible control u such that J(u) < oo and for the exponential 
stability of the closed-loop system we need the following two hypotheses. 

(HI). The system (3.1) is stabilizable in the sense that there is a feedback operator 

K E C(H, U) such that the closed loop semigroup Sk{1) E £(//) given by 

(3.5) S K {t)z = S(t)z + f S(t — s)BKSk{s)z ds 

Jo 

for all t > 0 and z (E H decays exponentially. 

(H2). The system (3.1) is detectable in the sense that there exists an operator F E 
C(Y,H) such that the output injection semigroup S F (t) E C(H) given by 

(3.6) S F {t)z = S(t)z + f S F (t - s)FCS(s)z ds 

Jo 

for all t > 0 and z E H decays exponentially. 

Remark 3.1. [28, pp. 134-135] (i) If (HI) is satisfied, then for any z 0 E H , there is an 

admissible control u 20 (-) E L 2 (0,oo',U) such that J(u zo ) < oo. 

(ii) Let (H2) be satisfied. Then for any z 0 E H and u(-) E T 2 (0, oo; U) with J(u) < oo, 
z(t) defined by equation (3.2) is in T 2 (0, oo; H). 

Now we state the following fundamental results for (LQR) problem, see [4,13]. 

Theorem 3.2. Let hypotheses (HI) and (H2) be satisfied. Then there is a unique optimal 
control u( ) E L 2 ( 0, oo; U ) for the linear quadratic regulator (LQR) problem and u(-) is given 
by the feedback law 

(3.7) u(t) = n z(t), t > 0 

where z(t) is the corresponding optimal trajectory and n E £(H) is the unique nonnegative 
self-adjoint solution of the algebraic Riccati equation 

(3.8) A t nz + UAz-nBR-'B*nz + C*Cz = 0 

for every z E 'B(A). Moreover, the closed-loop semigroup Sn(Q G tC,(H) generated by the 
operator A — BR~ l B*Ii decays exponentially. 
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Remark 3.3. The algebraic Riccati equation (3.8) can be defined for all 2 £ H , since the 
right hand side of equation 

(3.9) UAz = -A*nz + nBR- 1 B*Ilz-C*Cz 

is well-defined for all 2 £ H, T>(A) is densely embedded in H, and hence we can extend the 
left hand side of equation (3.9) to 2 £ H continuously. 

We now consider the control problem, and its linearization, for Burgers’ equation. Define 
the operator A e , as in Section 2, by A t <j> = t<j>" for all (j> £ D(A C ) = 0,^) ft Hq(0,£). For 

the control input operator B and the observation output operator C we consider the Hilbert 
spaces H = L 2 ( 0,£), U = R and Y = R fc . Assume that B £ C(U,H) and C £ C(H,Y) are 
defined by 

(3.10) Bu = b(-)u and Cz = ( 2 (^ 1 ), • • • , £(£*)), 
where fe(-) £ H, u £ U , and x, £ (0,^), 1 < i < k, are defined by 

| rxi+6 

(3.11) z(xi) = — z(x)dx. 

Zv J*i- 6 

In equation (3.11), 6 > 0 is chosen so that (x; — <5, x,- + <!>) C (0,£) for all 1 < i < k. 

Consider the following linear control problem. 

(LQR) a : Find u(*) 6 L 2 (0,oo]U) minimizing the weighted performance index 

J r°° 

1 {ll»Wlly + *IW*)l£}«’“<«. (o>0) 

0 

subject to the governing equations 

(3.13) j-z( t) = A ( z(t ) + Bu(t), 2 ( 0 ) = zq 

at 

(3.14) y(t) = Cz(t), t> 0, 

where A e , R, C are as above and R > 0. 

Remark 3.4. (i) Equation (3.13) corresponds to the system 
d d 2 

-z(t,x) = e-^z(t,x) + b(x)u(t), t > 0, 0 < x < i , 

u(<,0) = v(t,£) = 0, 
v(0,z) = u 0 (:r) . 
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(ii) For each i, 1 < i < k, z(x { ) given by equation (3.11) represents an average value of 
z(x) over a small neighborhood of X{. We can regard each x t , 1 < * <5 as the location of 
a “weight”. 

(iii) The weight function e 2at in the definition (3.12) of the cost functional J will play an 
important role in the exponential decay rate, see Theorems 3.9 and 3.10. However, it also 
gives rise to the question of existence of an admissible control u(-) such that J(u) < oo. 

For the control problem (LQR)^, we introduce an a-shifted control system [15]. Let 
z(t) = z(t)e Qi , u(t) = u(t)e Qt and y(t) = y(t)e at . We then have a modified linear control 

problem 

(LQR); : Find u G L 2 (0,oo;f/) minimizing the cost functional 

(3.16) J(ti)= / {||y(0lly + ^ll«( < )llir}^ 

Jo 

subject to 

4 z{t ) = (Ac + oil) z(t) + Bu(t), z(0) = Z 0 , 

(3.17) dt 

m = Cz{t ), t > 0. 

The solutions for the system (3.17) again are taken as mild solutions. If we solve the 
problem (LQR)„ and apply 

(3.18) u a (t) = e~ at ii{t) (t> 0) 

to the original control system (3.13)-(3.14), then the resulting optimal trajectory z a (t) will 
satisfy the inequality 

(3.19) \\z a (t)\\ H < Me~ at \\z 0 \\H, 

where M > 1 is a constant and a > 0 is the desired degree of stability. 

Remark 3.5. A discussion of the “a-shifted” problem for finite dimensional systems first 
appeared in [1]. Anderson and Moore showed that, for finite dimensional systems, the control 
problem (LQR) a is “equivalent” to (LQR); in the following senses : 

(i) The minimum value of J defined by equation (3.12) is the same as the minimum value 
of J given by equation (3.16). 
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(ii) If u(t) = g(z(t )) is the optimal control for (LQR); for some function g , then u(t) = 

e ~ at d(z(t)e at ) 

is the optimal control for (LQR) a and conversely. 

We now apply Theorem 3.2 to the problem (LQR); in order to obtain an optimal control 
u e L 2 (0, oo; U) for (LQR);. 

Remark 3.6. (i) From Remark 2.2, the spectrum cr{A t -f al ) of the infinitesimal generator 
A t +al consists of all eigenvalues A a>n , n = 1, 2, • • • , given by A a>n = a-en 2 n 2 /£ 2 and for each 
n, n = 1,2, •••, the eigenfunction 0 a>n corresponding to A Qi „ is given by <f> Q , n (x) = V^sm^fx. 

(ii) We are interested in the stabilization problem for the system (3.15) with small 
viscosity t = A. > o, i. e ., high Reynolds number. Let a 0 = If a > a 0 , then there is at 
least one positive eigenvalue A q ,i = a - Moreover, if a > a 0 and 0 < 7 < e, then the 
first eigenvalue o satisfies a — ^5- > a — > 0 and hence A ^ + al will 

have at least one positive eigenvalue for all 0 < 7 < e. This will become important when 
we apply feedback laws computed from low Reynolds number to investigate the closed-loop 
response of the nonlinear Burgers’ equation for high Reynolds numbers (see Example 3 in 
Section 4). 

Let a > a 0 = be given and let 

(3.20) n a = max { n G N : X a<n = a - —■* > 0 }. 

Since A c is self-adjoint (see Remark 2.2) and the set { <j>a,n • w = 1,2, •••} is a basis for 
H = L 2 ( 0,Q, we can identify z € H with the sequence {< z, </>«»,„ > }„ g n- Assume that 
b n G U and c n G Y satisfy 

OO 

(3.21) Bu = {< b n ,u >}„ € n and Cz = y^c n < z, <j) a>n > 

n= 1 

with zr=i Holler < 00 and ||c n | \y < 00, see [28, pp. 137-143]. 

The following lemma is an application of stabilizability and detectability results of 
Pritchard and Salamon [28, Section 4.2]. 

Lemma 3.7. For each n = 1,2, •• • , n a , let 

(3-22) „ = { — : i = 1 2, ■ - • , n — 1 } , 

n 
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Then the following statements hold. 

(a) b n =< b(-),<f> a , n >± 0 for all n = 1,2,- •• ,n a , if and only if the system (3.17) is 
stabilizable in H . 

(b) If 8 > 0 satisfies 6 < then for each n = 1,2, • • • ,n a , there exists at least one 
Xi , 1 < i < k, such that x { £ X a , n if and only if the system (3.17) is detectable through 

C G C(H,Y). 

Proof: (a) From Remark 3.6, we know that the spectrum a(A ( + al) of A ( + al consists 
of all eigenvalues A a , n = a - en 2 w 2 /£ 2 . Thus, for all n, n > n a + 1, we know that A a , n < 0. 
Let H u be the linear span of eigenfunctions <f> ail ,- ■ • , ^ t „„. Then the dimension of H u is n a 
and hence the system (3.17) is stabilizable if and only if the projection of (3.17) onto H u is 
controllable, or equivalently, if and only if b n =< b(-),(f> ain > ^ 0 for all n, n = 1,2, • • • , n a . 
For (b), let c n = (c nl , • • • , c nk ), n = 1,2,---, be defined by equation (3.21). Then we 

have 

\/2i nirxi . nirS 

(3.23) c ” = ^ Si ” T Sm T 

for 1 < i < k. By the dual statements of (a), the system (3.17) is detectable through 
C 6 C(H, Y) if and only if c n / 0 for all n = 1, 2, • • • , n a . Hence, (b) holds. 

Remark 3.8. If n a = 1, then X„,i is the empty set. 

We now return to the original control problem (LQR) a . The following theorem is the 
main result for our control problem (LQR) a . 

Theorem 3.9. Let a > a 0 be given. Suppose that b(-) G H = L 2 (0,t), 6 > 0, 1 < i < 

fc, satisfy the conditions (a) and (b) of Lemma 3.7. Then there is a unique optimal control 
UaA’) e F 2 ( 0,oo;R) for the problem (LQR) a such that 

(3.24) UaAt) — n a ,e ^a,t(f), t >0 , 

where z a , e (t) is the corresponding optimal trajectory and n a>£ 6 C(H) is the unique non- 
negative self-adjoint operator satisfying the algebraic Riccati equation 

(3.25) (A e + a/)*n at£ z + II a , £ (/l £ + al)z - I^SiT 1 B*Il ail z + C*Cz = 0 

for every z G T>(yl £ ) = H 2 { 0,0 n Moreover, the closed loop semigroup Sn o> «(0 e 

C(H) satisfies the following stability property 

(3.26) ||S„ a ,.(()IU W <Me-<” + “" 
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for some constants M = M(a , e) > 1 and u> — u(a, c) > 0. 

Proof: By Lemma 3.7, we know that the a-shifted control system (3.17) satisfies all 
hypotheses (HI) and (H2) (with z(t), y{t), u(t ), A , S(t ) and J replaced by z(t), y(t), u(t), 
At + al, S(t) and J, respectively ). Hence, by Theorem 3.2, there is a unique optimal control 
w(f) for (LQR) a and the corresponding closed-loop semigroup S(t) decays exponentially, i.e., 

(3-27) IlSWIk*) < Me'"' t> 0 

for some constants M = A/(a, e) > 1 and u> = tt>(a, e) > 0. Moreover, u(f) is given by 
(3.28) 6(0 = -RT'BrUaJit) 

where z(t) is the resulting optimal trajectory for the a-shifted system (3.17) and n a , e is 
the unique nonnegative self-adjoint solution of equation (3.25). Since the semigroup S^t') is 
generated by A t -\-otI SR 1 5*n Qje , the infinitesimal generator of the closed-loop semigroup 
S n a ,M for the original system (3.13) is A c - BR~ x B*Tl o t . Hence, 5n„ it (0 = S(t)e~ at and, 
by the relation (3.27), 5'n a< (0 satisfies the inequality (3.26) with M = M. Moreover, the 
optimal control u Q , e (<) for (LQR) q is given by the formula (3.24), since u Q , t (0 = it(t)e~ ai — 
—R l B B a ^z{t)e at = —R 1 B*Y[ Qyt z at (t), where z a ^(t) = z(t)e~ ai is the corresponding 
optimal trajectory for the original system (3.13). This completes the proof. 

The optimal control u a f (.) € T 2 (0, oo;R) obtained in Theorem 3.9 is given by the 
feedback law (3.24). Note that u Q ,t(f) depends on a and t. Define the feedback operator 
K a , t €C(H,U) by 

(3-29) K a>t = -R~ l B* n o , e . 

Then the optimal control u a>t (t), t > 0, is given by 

(3-30) *aA*) = K***aAi) 

and the infinitesimal generator for the closed-loop semigroup 5n at (f) is 

(3-31) A ( + BK ac = A t - BR-'B'H^. 

Recall that H = T 2 (0,^) and U = R. Thus, by Riesz representation theorem (see, e.g., 
[11, p. 13]), there is a unique feedback gain function k a>c (-) 6 L 2 (0,^) such that 

(3-32) I< aiC z= f k aiC (s)z(s)ds 
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for all z G L 2 (0,£). 


We now have the following main result of this section. 

Theorem 3.10. Let a > a 0 be given. Suppose that b(-) € H = L 2 (MK <5 > 0 and 
x u 1 < i < k, satisfy the conditions (a) and (b) of Lemma 3.7. Let k a>t {-) G H be the 
linear feedback gain function defined by the formula (3.32). Then there exist constants 
p = p(a,e) > 0 and M = M(a,e) > 1 such that for any initial data z 0 (-) € #o(M)> Wlth 
||z 0 || H i < 2 ^, the controlled Burgers’ equation 

z(t,x ) = e -^—z(t,x) - z(t,x)—z(t,x) + b(x) / k a<c (s)z(t,s ) ds 
at ox 1 ox jo 

(3 ' 33) z{t,0) = z(t,e) = Q, 

z(0, x) = z 0 (x) G Hq(0,£) 

ha s a unique (strong) solution and the solution *(*)(•) = z{t, •) satisfies the following stability 
property 

(3.34) IPWIIffo 1 < 2Me~ at IMOIIhJ- 

Proof: Let the operators A„ B, C and I< a , t be given by equations (2.4), (3.10) and 
(3.29). Define the nonlinear function / : H^( 0,£) — ♦ L 2 ( 0,f) by 

(3.35) f( z ) = BK °,t z ~ zz ' 

where ' = — . Then, the map / satisfies the hypothesis (F) in Section 2, since for any 

dx ' 

z i^GtfoW), 

( 3 . 36 ) ||/(2i) - f(z?) ||l»(o,/) < ( \\BK a , t \\ C (H) + II^iIIh* + INI h\ ) \\ z i ~ z 2||//i- 

Note that the operator BI< a , ( is bounded on the state space H = L 2 {0,1). Thus, by Theorem 
2.1, we have a unique local (strong) solution of equation (3.33). 

Let Sx a <(*)) t > 0, be the analytic semigroup on H generated by the operator A t + 
BI< a<e . Then, by Theorem 3.9, S Ka( {t) satisfies the inequality 

||^ a , e (02||//<M a e-(° +w)t ||z|| // 


(3.37) 
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for all z € H and for some constants M a = M a (a, c) > 1 and u> = w ( q , c ) > 0. Let 
a < P < a -|- <jo. Then, there is a constant M a = M a (a,e,fl) > 1 such that 

(3-38) I|Sk„«*II»J 

(3-39) l|S*„W*llflS 

for all 2 € Hq. Let M = max{ M a , M a } and choose p > 0 with 0 < p < Then it is 

easy to see that 

° 1 1 

(3.40) pM — e -W- a > ds < 

Jo v s 2 

Thus, by arguments similar to those in the proof of Theorem 2.7 together with inequalities 
(3.38)-(3.39) and the expression 

( 3 - 41 ) *(0 = S K a ,'(t) z o + f S Ka t (t - s)g(z(s))ds, 

JO 

the unique global solution 2(t)(-) = z(t,-) for the controlled Burgers’ equation (3.33) exists 
and satisfies the inequality (3.34), where g{z(t)) = —z(t)z'(t). 

4. Approximation and numerical results 

In Section 3, we considered a linear quadratic regulator problem (LQR) a to obtain a 
desired degree of stability for the solution of the closed-loop Burgers’ equation. In this 
section we consider an approximation scheme for (LQR) a and give some numerical results. 

We first introduce an abstract approximation scheme for the problem (LQR) based on 
the results of Banks, Kunish [4] and Ito [19] and then apply the scheme to get the optimal 
control «<»,((•) € T 2 (0, oo; R). Throughout this section, we assume that R = /, the identity 
operator on the control space U, (see Section 3). For more approximation schemes for the 
linear regulator problems, see [2,4,14] and references given there. 

We formulate a sequence of approximate regulator problems and present a convergence 
result for the corresponding Riccati operators. Throughout this section, we use superscript 
N in the designation of subspaces, operators and matrices in the A-th approximating system 
and control problem, like H N ,A N , B N , etc. Hence the superscript N indicates the order of 
approximation. 
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Let , N — 1,2, • • • , be a sequence of finite dimensional linear subspaces of H and 

pN . fj fjN t ^ e canonical orthogonal projections. Assume that S N (t) is a sequence 

of Co-semigroups on H N with infinitesimal generators A N G C(H N ). We then consider the 
family of regulator problems: 

(LQR) n : Minimize J N (zg,u) over u G L 2 (0,oo;C) subject to the control system 

(4.1) z N (t) = S N (t)zg + J S N {t-s)B N u(s)ds, 

(4.2) y N (t) = C N z N (t ), 

where z N ( 0) = z A = P N Zq and 

roc 

(4.3) J N (zo, u )= / { ll^COIIr + IKOIlM dt - 

Jo 

Remark 4.1. If, for each N, [A N ,B N ) is stabilizable and ( A N ,C N ) is detectable, then, 
by Theorem 3.2, there is a unique optimal control u N (t) for the finite dimensional problem 
(LQR) n and it is given by 

( 4 . 4 ) u N (t) = -{B"yn s sX(t)z?, 

where Sj({t) is the C 0 -semigroup on H N generated by A N - B N (B N )*U N and ll‘ v G C(H N ) 
is the unique nonnegative self-adjoint solution of 

(4.5) (A N )-n K + n"A w - n w B' v (B'')-n'' + (c N yc N = o. 

For the finite dimensional approximation systems, it is not clear that ( A ,B ) is sta- 
bilizable even if the original system (A, B) is stabilizable. Similarly, it is not clear that 
the detectability property of (A, C) is preserved under the finite dimensional projections. 
Another question we have to consider is the convergence of approximates U N and u A (t) to 
the infinite dimensional solutions II and u(t), respectively. For these reasons, we need the 
following assumptions. 

Let S N (t) = e ANt , t > 0. 

(Al): For each z G H, S N (t)P N z — - S{t)z and S N {tyP N z — ► S{t)*z, where the 
convergences are uniform in t on bounded subsets of [0, oo). 

(A2): (i) For each u G U, B N u — + Bu and for each z € H, ( B N )*P N z » B’z. 


19 



(ii) For each z G //, C N P N z — ► Cz and for each (C N )*y — » C*y. 

(y!3): (i) The family of systems (A N ,B N ) is uniformly stabilizable , i.e., there 

exists a sequence of operators K N G such that sup ||/T^|| < oo and 

< M.e-",', <>0, 

for some positive constants Mi > 1 and u)\ > 0 which are independent of N. 

(ii) The family of pairs (A N ,C N ) is uniformly detectable , i.e., there exists 
a sequence of operators F N G £(K, //^) such that sup ||F lAr || < oo and 

|| e (X» + F»C»X^ ||£(Jf) < ( > 0, 

for some constants M 2 > 1 and W 2 > 0 which are independent of TV. 

Remark 4.2. (i) The condition (y43)(ii) is a relaxation of the coercivity assumption in [4] 
(see also [19, p.3]). 

(ii) Suppose that B N = P N B and C N — CP N . Then (.42) holds, since it follows from 
(41) that P N z — » z for all z £ H. 

By simple modification of results from [19, Theorem 2.1] and [4, Theorem 2.2], we have 
the following fundamental convergence results. 

Theorem 4.3. Let (A, B) be stabilizable and (4,C) be detectable. Suppose that (41) — 
(43) are satisfied. Then, for each N, the finite dimensional algebraic Riccati Equation 

(4.5) admits a uique nonnegative self-adjoint solution 11^ such that sup{ ||n^|[£/^) : N = 
1, 2, • • • } <00 and 

(4.6) — * n z 

for every z 6 H. Moreover, there exist positive constants M 3 > 1 and u> 3 (independent of 
N) such that 

(4.7) | |e (A»- e »(B''|-n»),piV|| £(/f| < ( > 0 


For the uniform stabilizability and detectability assumptions (43) we introduce a sesquilin- 
ear form a t (',') : V x V — * C defined by 


(4.8) 



ez'(x)w'(x)dx , 


z,w 6 V , 
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where V = Hq(0,1). Note that, to allow the use of the theory of sectorial operators and 
sesquilinear forms in discussing the spectra of various operators, we assume in defining a t (-, •) 
that the functions in V are complex valued. It is easy to see that the sesquilinear form a £ (-, •) 
is P-coercive [36, p. 274], i.e., 

(4.9) |a £ (.z,u>)| < e||z||v |Mlv> (continuity), 

(4.10) 'Rea t {z,z)-\-'i\\z\\ 2 H >€.\\z\\y, (Garding’s inequality) , 

for all z,w £ V and 7 > e > 0. Furthermore, it follows from the bounds (4.9) and (4.10) 
that there exists, in a unique manner, an operator A c € C(V, V*) such that 

(4.11) a ( (z, w) =< — A ( z, w >v,v and a ( (z,w) =< —A*w, z >v,v 

for all z,w £V (see, e.g., [36, pp. 271-275]). 

Turning next to specific approximations for ( LQR) a , we divide the unit interval [0,1] into 
N + 1 equal subintervals to get [27, ar,-+i], £,• = 77+1 ’ * = 0, 1, ■ • • ,N. For each i, 1 < * < N, 
let h^(x) denote the linear spline basis function defined by 

{ (N + l)(x - Xi-i), £i-i<£<£,- 
~(N + l)(x - x, + i), £,<£<£,+ 1 

0 , otherwise. 

Let H N be the ^-dimensional finite element space given by 

N 

(4.13) H N = {'$2 zih?(x):zi eft, i = 1,2,- ,N}. 

t = 1 

Then we have a sequence of finite dimensional (real) subspaces H N C V, n = 1 , 2 , •••. 
Moreover, it is well-known [31], [19, p. 15] that the family of H N satisfies the following 
approximation condition: 

(APP) : For each z £ V, there exists an element z N € H N such that 
\\z — z N \\v < e(N), where e(N) — > 00 . 

Let P : H —y H N be the canonical orthogonal projection onto H N . Then, from the 
approximation property (APP), it is a trivial matter to see that P N z -> 2 as JV -+ 00 , for z 6 
H. For the finite demensional regulator problem ( LQR )% we choose 

(4.14) B n = P n B and C N = CP N . 
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Then the conditions (A2)(i),(ii) follow from Remark 4.2. 

To obtain a representation A ^ of A ( on H N , consider the restriction of the sesquilinear 
form a t (-, •) to H N x H N . We then have a representation A ^ of A t satisfying 

(4.15) a c (z,w) =< —A^ z,w > and a t (z,w) =< — (A^)*w,z > 

for all z,w G H N . Equation (4.15) follows from the fact that H N is a real Hilbert space. We 
know also that A ^ = (A^)*, since A e = A*. 

Remark 4.4. Since H N C H, by equation (4.15), it is easy to see that for any A G < 7 (A^), 
He\ < —7 < — e. 

Let S N (t) be the Co-semigroup generated by A^ . Then the conditions (Al)(i),(ii) fol- 
low from the results of Banks and Kunish [4, Lemma 3.2]. Note that S N (t) = (S N (t))*. 
For the condition (A3)(i) we need a certain preservation of exponential stabilizability under 
approximation ( (POES) in [4]). The following result is taken from [4, Lemma 3.3]. 


Theorem 4.5. Let (A e , B) be (exponentially) stabilizable. Suppose that the approximation 
condition (APP) holds. Then the approximations defined through equations (4.14)-(4.15) 
satisfy the condition (A3)(i), i.e., the family of pairs (A?,B N ) is uniformly stabilizable. 


By the dual arguments of Theorem 4.5 we can see that the condition (A3) (ii) holds 
under the assumption that (A t ,C) is detectable. We summarize our discussion up to this 
point as the following theorem. 

Theorem 4.6. Let (A t + al , B) be stabilizable and (A c + aI,C) be detectable. Let A^, 
B n , C N be defined as in equations (4.14) and (4.15). Then we have 

(4.16) n£ £ P N 2 — ► U a<( z, zeH , 

and 

(4.17) S N (t)P N z — + S{t)z , zeff, 


where the convergence is uniform in t on bounded subsets of [0, 00 ), P N is the orthogonal 
projection onto H N , and n^ £ satisfies 


(4.18) (A? + a/ N )*n^ + n"(A" + cd N ) - n^B N (B N yn^ + (C N )*C N = 0 . 
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Remark 4.7. Note that S(t) = S(<K“‘, = S N {t)e~ at , where S(t), S N {t ), S(t) and 

S N {t) are semigroups generated by >4 £ , A?, A e + al and A? + aI N , respectively. 

Next, consider the matrix representations of operators on the space H N . Let the ap- 
proximate solution z N (t,x) of z(t,x) on H N be given by 

(4.19) z N (t,x) = £ zf'Wl' (*) 

1 = 1 

for some zf(t) 6 R, i = 1, * • • , N. Then, from equations (4.14) and (4.15), we have a finite 
dimensional ODE system 


(4.20) 

where {z N (t)} = [z^ ( t ), • • • ,^(0] T > 

[< G N } = [< >] n *n 


lG N ]±h N (t)} = [A?]{z N (t)} + {B N Mi), 

at 


(4.21) 


6(N + 1 ) 


4 10 0 

14 1 0 

0 14 1 


(4.22) 


K] = e(N + 1) 


0 
0 

-2100 
1-210 
0 1-21 


0 

0 


... 0 

... 0 

... 0 


1 4 1 

0 1 4 


J NxN 

0 

0 

0 


1 -2 1 

0 1 -2 


NxN 


(4.23) {B N } = [< 6, h? >, < b, h” >, • • • , < ^ h% >] T , 

where < b, hf >= £ b(x)h?{x)dx, 1 < j < N. Since [G N ] is invertible, by multiplying 
[G"]” 1 to both sides of (4.20), we get 


(4.24) 
where 

(4.25) 


~{z N (t)) = «l(z N (t)} + {£")«(<), {*"(0)} = {*"}. 

dt 


| A N 1 = [ G N ]- , [^] 1 { b ») = 
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and {z*} = [<?*] '[< z 0 ,h? >,•••,< z 0 , h% >] T . 

Next, consider a representation C N of the operator C on H N . It is easy to see that 
C N : H n — > R k is given by 


( 4 - 26 ) [C N ] = [hf (5.-)W. 

where h 1 - (x.) = ^ A^(x)dx, 1 < t < fc, 1 < j < N. 

Finally, we have a finite dimensional Riccati equation (4.18) and the corresponding 
feedback gain operator K^ t given by 


(4.27) K” f = -(B")-n« . 

Therfore, the closed loop system (4.24) can be represented by 


(4.28) 


-{z"(t)} = (A? + B N Kl){z N (t)}, {z»} = {zg}. 


Now, we discuss an algorithm for finding the unique nonnegative self-adjoint Riccati 
solution for equation (4.18). We employ the Potter’s method [30] to obtain II J e . The first 
step in Potter’s method is to form 2 N x 2 N matrix 


(4.29) 


M n 


(a? + ai N y (c N yc N 

B n (B n )* ~(A? + aI N )/ 


Next, find all eigenvalues and eigenvectors of M N and form the matrix 


(4.30) 


Z N 


Q f 
Q!\ 




where the columns of Z N are the eigenvectors of M N corresponding to the eigenvalues 
with positive real part. When eigenvalues occur in complex conjugate pairs, so do the 
eigenvectors. In this case, the real and imaginary part of the eigenvector each forms a 
column of Z N . Finally, the solution to the Riccati equation (4.18) is given by the formula 

n" = Q?(Q?r l . 


Remark 4.8. From the numerical results we found that the Riccati solution operators n^ f 
blow up when the viscosity e > 0 goes to 0 for fixed a > 0. Also, when a goes to infinity 
with e fixed the same phenomenon has been observed. 
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Finally, the finite dimensional approximation for the controlled nonlinear Burgers’ equa- 
tion (3.33) is given by 


(4.3i) jyw = «-»"] + {fl'hraxAo} + /"asm. <*»> = 


where [/l( v ], ( B N ) are defined as in equation (4.25), and 

/"({^(f)}) = [G K ]-'/ N ({^(i)}), 

^( t )^(0 + (^ W ) 2 


/"({a* (()» = - 


-( 2 «( t )) 2 - (*?(<))(*?«) + (' 4 '( t ))( 4 < t )) + (^(()) 2 


where [G"] is defined by equation (4.21). To solve the nonlinear ODE system (4.31) we use 
the 4-th order Runge-Kutta method [33]. 

Throughout the rest of this section, we discuss how our results work for relaxation of 
“steep” gradient of the solution for Burgers’ equation through numerical experiments. 

For numerical examples, the length £ for space domain, the Reynolds nummber, Re, the 
initial function z 0 ( ■) € H^(0, 1) and the control input function &(•) € T 2 (0, 1) will be chosen 
as 1, 60, sin nx and e x , respectively. Thus the governing equation is given by 


(4.32) 


!*(«.*) = k - z{t ’ x) l z(t ' x) + e * i f WsWM) ds 


*(t,o) = z(M) = o 

z(0, x) = sin 7 rx, 


where the feedback gain function k a<e {-) € L 2 { 0, 1) will be determined by the desired degree 
a > 0 of stability and the action of output operator C. 

The “robustness” of the feedback controller exhibited, for example, in Figure 2, will 
be discussed in Example 3. For this particular example, Reynolds numbers 60, 80, 100 and 
120 are chosen. 


Remark 4.9. (i) From the numerical experiments, we found that if Reynolds number, 

Re, is less than 60, then the diffusion phenomena dominate convection phenomena. In this 
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case, the formation of steep gradient due to convection term — z(t, x)-^z(t, x) of the open- 
loop solution, i.e., k Q<t (-) = 0 in equation (4.32), is not clear. But, for Reynolds number 
greater than 60, the open-loop solution creates “sharp” gradient in finite time, (see Figures 
1, 10, 12 and 14). Of course, the solution dies out eventually, because of the diffusion term 

i x ) • 

(ii) The control input function b(x) = e x is defined for all x G [0,1]. Thus, the 

feedback control acts on the whole domain [0,1]. But, one can choose any L 2 -function 
K‘) G T 2 (0, 1) satisfying the stabilizability condition in Lemma 3.7. In fact, b(x) = e x 
satisfies the stabilizability condition for any desired degree of stability a > 0, since the 
coefficients b n , n = 1,2,--*, representing input function &(•) are not zero, i.e., b n = < 
6(-),sin nx >l 2 (o,i)— e x sin7 rxdx ± 0 for all n = 1,2, •••, (see Lemma 3.7). 

(iii) The initial function Zq(x) = sin7rx is chosen for our numerical experiments. Other 
typical //^-functions such as the “hat function” defined by 


(4.33) 


z 0 (x) 


2x, x G [0, |] 

-2x + 2, x G [|, 1] 


can be used for initial data. But, we found that the solution of Burgers’ equation (4.32) with 
initial data 2o( a; ) replaced by the hat function has almost similar phenomena, such as the 
creation or relaxation of steep gradients, as those of sloution with initial data z 0 (x ) = sin nx. 

To show trajectories of open- loop and closed- loop solutions, the order N of approximation 
is chosen as N = 32 for both cases. And the corresponding trajectories from time t = 0.0 to 
t = 1.0 will be shown. The convergence of the feedback gain functions k a<c (-) G L 2 (0, 1) will 
be shown for N = 8, 16, 32, 64 and 128. 


Example 1 . (Bounded Input/Output) 

The observation operator C G £(L 2 ( 0, 1),R 3 ) for this example is given by 

(4-34) C(z) = (5(0.3), 5(0.5), 5(0.75)), 


where z(x) is the average value of z(-) G L 2 (0, 1) in a small neighborhood of x, x = 
0.3, 0.5, 0.75, and defined by equation (3.11), 5(x) = ^ z(s)ds. Here, 6 > 0 is chosen 
so small that each open interval (x - S,x + 6) is contained in the whole domain (0, 1). The 
desired degree a of stability are chosen 0.3 and 0.6 for Figures 3 and 5, respectively. For 
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both cases, n a = max{ n £ N : a — ^ n 2 n 2 > 0 } — 1 and hence the set X at i defined in 
Lemma 3.7 is empty. Thus, all assumptions in Theorem 3.9 are satisfied. 

The feedback gain functions k a<c (-) are given in Figures 2 and 4. From these plots, it is 
easy to see that control action is concentrated on the location of sensors. This phenomenon 
is natural, since the optimal control is obtained to minimize the cost functional J defined by 
equation (3.16) whose first term ||j/(0llv' = ||C2(0 IIr 3 = X);=i l^(^*)| 2 > where Xi — 0.3, 0.5 
and 0.75 for i = 1,2 and 3, respectively. The corresponding closed-loop trajectories are 
shown in Figure 3 (for a = 0.3) and Figure 5 (for a = 0.6). From Figures 3 and 5, we can 
see how the controllers contribute to stabilization of the steep gradient as well as the solution 
itself. 

Example 2. (Identity Output Operator) 

For this example, we take the identity operator / on Z/ 2 ( 0, 1) for the output operator C. In 
this case, the ouput space Y is L 2 ( 0, 1). The convergence of gain functions and corresponding 
closed-loop trajectories for a = 0.3 and 0.6 are shown in Figures 6-7 and 8-9, respectively. 
Since the observation operator is identity, this example gives the information about the 
maximal control action. We note the following observation concerning the convergence rate 
of gain function. Theoretically, the rate is O(^) [19, p. 15]. But, in this example, the rate 
seems to be faster than 0(£), (see Figures 2 and 4). Another observation is concerned with 
the location of maximal control action. The location moves to the left portion of domain as 
the degree of stability a > 0 increases. In other words, we should put more action on the 
front part of domain to get a higher exponential decay rate a > 0. (See also Figures 2 and 

4 ). 

Example 3. (Robustness) 

In this example we show the robustness of the feedback controller showed in Figure 2. 
The feedback controller is obtained from the control system with Re=60, a = 0.3, b(x) = t x 
and the ouput operator C defined by equation (4.34). 

Figures 2, 10, 12 and 14 show open-loop trajectories for Reynolds numbers 60, 80, 
100 and 120, respectively. The corresponding closed-loop trajectories are shown in Figures 
2, 11, 13 and 15, respectively. The order N of approximation is chosen as N = 32 for 
Re=60, 80, 100 and N = 64 for Re=120. From these examples, it is easy to see that 
the feedback controller obtained for Re=60 stabilizes the steep gradient of the solution for 
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Burgers’ equation with various Reynolds numbers. However, we see that the sharp gradient 
is relaxed slowly as Reynolds number increases, (see Figures 3, 11, 13 and 15). Note that, 
from Theorem 3.10, the closed-loop solution z(t) of the nonlinear system (3.33) satisfies the 
stability property 

(4.35) IMOIIfl’ < 2Me-"' |M-)||„, . 

Although the exponential decay rate a is independent of Reynolds number, the constant 
M = M(a, Re) depends strongly on the Reynolds number, Re. 

5. Conclusion 

In this paper, we considered a feedback control problem for a nonlinear equation, in 
particular, Burgers’ equation. The method consists of linearization of the nonlinear equation. 
We used the linear quadratic regulator (LQR) problem to find optimal feedback gains. The 
linearized equation is the heat equation. It was also proved that, under appropriate selection 
of the input functions and “weights location”, the LQR problem for the linearized problem is 
detectable and stabilizable. We then analysed a “shifted quadratic cost” to construct gains 
which produce a fixed decay rate. In particular, we showed that the closed-loop system 
satisfies the inequality 

\\z(t y 0]z o )\\ H i < M(t)e Ikolltfi, 

where a > 0 does not depend on the Reynolds number, but M(e) does, (see Theorem 3.10). 

We also developed a numerical scheme for computing the feedback functional gains. 
Several numerical experiments were performed and the following observations were made: 

1) The functional gains depend strongly on the “weight location”. For example, if the 
output operator C is given by 

C{z) = (5(.3),z(.5),z(.75)), 

then the gain function is concentrated on the weight location, 0.3, 0.5 and 0.75, (Figures 2 
and 4). 

2) The closed-loop nonlinear system is stabilized (as predicted) by linear feedback laws. 
Moreover, the steep gradients (for e « 0) are smoothed out by feedback. 

3) To test the “robustness” of the feedback control law, one experiment was performed. 
We obtained the functional gain fc^ e (-) from the control system at the Reynolds number, 
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Re=60, and applied it to the closed-loop system 


d , N 1 
( 5 - 1 ) di Z ^' X ^ ~ 


d 2 


Re&r 2 


r\ rl 

z(t,x) - z(t, x )^ z (^ x ) + e " Jf k lA s ) z (^ ^ ds 


at Re=80, 100 and 120. The closed-loop responses are shown in Figures 11, 13 and 15. 
Although the performance was decreased, the system (5.1) was still stabilized and smoothed 
out. These results provided some insight into the possibility of using linear feedback laws 
for nonlinear distributed parameter systems. 
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Figure 11 Closed Loop (R« = 80, feedback by k^( J) 


Figure 12 Open Loop (N -32, Re=100, *o(*) - »in»*) 
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Figure 14 Open Loop (N-M, Re 120 zo(r) = mwt) 
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